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ABSTRACT 


In  this  special  technical  report  we  present  fundamental  studies  on 
Instability  analysis.  When  an  electromagnetic  wave  linearly  Interacts 
with  an  active  medium  or  system,  there  can  occur  basically  two  types  of 
Instabilities.  One  of  them,  In  which  the  wave  grows  Indefinitely  with 
time.  Is  called  absolute  Instability,  and  the  other.  In  which  the  wave 
grows  spatially  as  It  travels.  Is  called  convective  Instability.  By 
conformal  mapping  of  the  dispersion  equation  from  cu  onto  the  k  plane  or 
vice  versa,  we  have  made  a  detailed  analysis  of  these  two  types  of  Insta¬ 
bilities.  General  procedures  for  analyzing  Instabilities  are  summarized. 
Effects  of  driving  sources  and  the  possible  occurrence  of  branch  cuts,  In 
the  k  plane,  of  the  dispersion  equation  are  also  discussed. 
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INTRODUCTION 


In  recent  years  the  geophysical  scientific  conununlty  has  shown  con¬ 
siderable  Interest  In  the  Idea  of  magnetospheric  plasma-lnJectlon  experi¬ 
ments.  It  Is  hoped  that  these  experiments  will  aid  in  the  study  of  the 
following  Important  problems: 

(1)  Particle  precipitation  and  the  physics  of  the  lower  ionosphere. 

(2)  The  details  of  wave-particle  interactions,  which  lead  to  the 
full  understanding  of  the  mechanisms  causing  natural  emissions 
and  the  control  of  artificially  triggered  emissions. 

(3)  The  physics  of  aurorae. 

(4)  The  feasibility  of  using  an  appropriate  controllable  amplified 
wave  for  space-to-space  and/or  space-to-ground  communications. 

In  all  the  above  problem  areas,  it  is  of  primary  importance  to 
understand  the  details  of  wave-partlcle/medium  interactions  and  the 
related  instabilities. 

When  a  wave  interacts  with  a  linearized  system,  two  fundamental  types 
of  instability  (absolute  and  convective)  can  occur.  The  detailed  distinc¬ 
tions  between  the  two  are  given  in  Section  III.  Due  to  various  applica¬ 
tions  of  instability  theories  in  both  the  space  experiments  and  the 
laboratory  plasma  devices,  such  as  fusion  and  beam-plasma  devices,  the 
subject  of  wave-plasma  instability  has  been  investigated  by  a  great  num¬ 
ber  of  workers  (listed  in  the  Bibliography).  The  instability  analysis  and 
criteria  and  the  computation  for  the  various  growth  rates  have  been  fairly 
well  developed.  However,  in  most  of  the  earlier  articles  on  the  subject, 
the  authors  give  only  brief  analyses  on  instability  theory.  Also,  the 


theorems  and  the  analyses  for  determining  the  nature  of  Instabilities  are 
scattered  in  the  various  Journal  articles  and  reports.  It  Is  the  purpose 
of  this  report  to  present  an  Integrated,  detailed  study  on  the  subject 
of  instability  analysis.  Starting  from  the  very  elementary  analysis,  we 
proceed  step  by  step  toward  the  conclusions  on  the  theory  of  Instabilities 

It  is  hoped  that  the  present  report  can  be  of  use  to  part  of  the 
geophysical  community  in  understanding  the  subject  of  wave-plasma 
(particle)  instabilities  and  their  related  applications  in  the  space 
plasmas. 

We  do  not  cite  the  references  In  the  various  sections.  The  impor¬ 
tant  publications  used  in  the  preparation  of  the  present  report  are 
listed  In  the  Bibliography. 


I I  DISPERSION  RELATION 


When  we  consider  the  problem  of  electromagnetic  radiation  or  wave 
propagation  in  a  uniform,  time-invariant  medium,  we  can  apply  the  Fourier 
transform  to  both  time  and  spatial  variables.  Using  linearized  source- 
free  Maxwell's  equations  combined  with  the  appropriate  dynamic  relations 
to  determine  the  constituent  parameters  such  as  permittivity,  c,  permea¬ 
bility,  ^l,  and  conductivity,  a,  we  can  derive  the  so-called  dispersion 
equation  D(cu,k)  =  0.  Wliat  we  customarily  know  of  the  dispersion  equation 
is  the  following: 

(1)  The  solutions  of  D(u),k)  =  0  determine  the  characteristic  wave 
modes  of  the  system.  When  both  cu  and  ic  in  these  solutions  are 
real,  the  characteristic  modes  are  then  representable  by  steady- 
state  plane  waves,  with  their  phase  varying  as  exp  ±  j(a)t  -  it  •  x). 

(2)  In  a  linear  regime,  when  we  perturb  the  system  by  a  small  source, 
the  response  r(x,t)  can  be  expressed  by  the  superposition  of  the 
characteristic  wave  modes. 

(3)  Whenever  D(u},k)  is  an  algebraic  polynomial  of  finite  order,  the 
system  can  support  only  a  finite  number  of  characteristic  modes. 
On  the  other  hand,  when  D(u),k)  is  a  transcendental  function  or 
an  Integral  form  (e.g.,  Landau  waves,  cyclotron  harmonic  waves, 
Bernstein  modes  in  warm  plasma)  the  system  may  support  an  infi¬ 
nite  number  of  characteristic  modes.  The  solutions  uj  (it)— here, 
both  03^  and  k  are  assumed  to  be  real— are  classified  in  terms  of 
the  various  branches  in  the  dispersion  diagram  (cu  vs.  it  plots). 

An  essential  feature  to  note  is  that  the  modes  uj  (it) ,  for  the 

n 

same  n  and  different  it’s,  as  well  as  for  different  n's,  are 
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linearly  independent.  Alternatively,  we  can  describe  the  wave 

types  0)  (Tc)  by  the  solutions  V.  (cb) .  For  the  case  of  both  co  and 
*1  n 

li,  considered  here  to  be  real,  there  is  a  one-to-one  correspon¬ 
dence  between  the  branches  o)  (S)  and  Hi  (co) . 

n  n 

(4)  The  phase  velocity  V  =  (a)/k)k  and  the  group  velocity  v  =  V-co 

P  B  k 

for  the  characteristic  waves  can  be  directly  derived  from  the 

dispersion  equation. 

The  above  understanding  of  the  dispersion  equation  is  sufficient  for 
analyzing  the  waves  in  simple  systems  Involving  passive  media.  Since  in 
these  simple  systems  it  is  not  possible  for  the  waves  to  draw  the  energy 
from  the  medium  or  the  system,  the  situation  for  the  wave  system  is  stable. 
In  a  stable  system,  the  characteristic  modes  determined  from  the  dispersion 
equation  are  either  steady-static  propagating  plane  waves  (both  cd  and  k 
are  real)  or  attenuated  plane  waves  (o)  real,  k  complex). 

When  we  deal  with  electromagnetic  problems  in  a  complex  system 
involving  active  media  such  as  warm  plasmas  and  energetic  beams,  it  is 
important  for  us  to  know  the  further  implications  as  well  as  mathematical 
properties  of  the  dispersion  equation. 

In  the  abovementioned  complex  system  it  is  possible  for  the  waves  to 
draw  energy  from  the  medium  or  system,  causing  the  wave  perturbation  to 
grow.  The  wave  system  is  thus  referred  to  as  an  unstable  system.  In 
analyzing  the  wave  properties  in  an  unstable  system,  several  additional 
implications  of  the  dispersion  equation  D(u),lS)  =  0  should  be  noted: 

(1)  D(a),k)  is  regarded  as  a  function  of  two  complex  variables,  o),  S. 

It  should  be  noted  that  in  general  there  is  no  one-to-one  corre¬ 
spondence  between  the  branches  fi  (oi  real)  and  co  (S  real).  How- 

n  n 

ever,  a  connection  can  be  made  through  a  process  of  conformal 
mapping  in  which  o)  and  li  are  regarded  as  complex  variables. 
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(2)  For  an  active  system,  there  arises  a  difficulty  In  Interpreting 
the  roots  with  it  complex  derived  from  solving  D(u),it)  =  0 
directly.  This  Is  because  we  cannot  distinguish,  purely  from 
the  sign  of  Im  k  In  exp  ±  jvcot  -k  •  x) ,  whether  the  wave  Is 
amplifying  or  attenuating,  since  we  do  not  know  whether  to  con¬ 
sider  Increasing  or  decreasing  values  In  spatial  coordinates. 

For  stable  systems,  the  complex  k  roots  cause  no  ambiguity 
because  on  physical  grounds  all  such  waves  must  be  attenuating. 
In  view  of  the  foregoing  arguments,  we  emphasize  that.  In 
unstable  systems,  solving  for  the  roots  with  ^  complex  and 
specified  cu  fran  D(cu,k)  =  0  has  little  significance. 

(3)  The  roots  of  D(a3,lc)  =  0  with  Im  cd  ^  0,  where  ">"  and  "<"  apply 
to  exp  -  j(cut  -Ic  •  x)  and  exp  +  j(a>t  -Ic  •  x),  respectively,  Indi¬ 
cate  that  the  active  system  Is  unstable.  Thus,  Instability  In 
some  sense  does  exist  In  the  waves.  An  ambiguity  does  not  arise 
In  the  complex  co  roots  here,  as  discussed  In  Itfcirn  2,  because  we 
always  consider  the  Increase  In  time  (causality). 

In  this  section  we  have  described  several  basic  Implications  of  a 
dispersion  equation.  It  should  be  mentioned  that  many  more  Implications 
than  we  have  listed  above  can  be  drawn  from  the  dispersion  after  a  detailed 
analysis.  In  fact,  a  sufficient  analysis  of  D(a),^)  =  0  can  produce  all  the 
Information  and  the  physical  quantities  regarding  wave  propagation  In 
either  a  stable  or  an  unstable  system. 

The  main  purpose  of  the  present  report  Is  to  provide  a  detailed  study 
on  the  problem  of  wave  radiation/propagation  In  an  unstable  system. 
Specifically,  we  examine  the  classification  of  Instability,  Its  mathe¬ 
matical  criteria,  and  the  computation  of  growth  rates.  All  these  entitles 
are.  In  fact,  contained  In  the  dispersion  equation.  One  method  of  recog¬ 
nizing  them  explicitly  Is  to  use  the  dispersion  equation,  causality  or  an 
appropriate  Initial  condition,  and  the  boundary  conditions,  and  to  perform 
a  detailed  Jtiialysis  on  both  the  complex  oj  plane  and  the  complex  Ic  plane. 
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Ill  CLASSIFICATICm  OF  INSTABILITIES 


Basically,  two  distinct  types  of  instability  can  occur  when  a  wave 
Interacts  with  an  unstable  system.  The  two  types  are  termed  absolute 
instability  and  convective  instability.  The  distinction  between  them 
becomes  apparent  only  after  a  spectrum  of  waves  is  superimposed  by  in¬ 
verting  the  Fourier  and  Laplace  transformations  (details  of  which  are 
given  In  the  next  section). 

We  consider  an  initially  quiescent  and  uniform,  time-invariant 
medium,  which  permits  us  to  look  for  solutions  such  as  exp  J(u)t-Ic  ‘  x) 
of  the  linearized  wave  equation.  The  frequency  and  wavenumber  are  con¬ 
nected  through  the  dispersion  equation  D(a),Ic)  =0.  If  Ic  is  real  and 
Im  u)  <  0,  this  solution  implies  that  the  amplitude  of  the  perturbation 
will  Increase  exponentially  with  time  for  any  x.  Hence,  within  the  limit 
of  the  linearized  theory,  the  system  is  unstable  and  no  steady-state  con¬ 
ditions  appear  possible.  To  examine  the  nature  of  increasing  amplitude 
of  the  perturbation  we  note  that  the  response  r(t,x)  can  be  expressed  by 

r(t,x)  =  (2rt)  ^  ^ doj  ^dS  i?(u),k)e^ ^  (1) 

L 

where  L  denotes  the  Laplace  path,  and  X(u),1c),  inversely  proportional  to 
D(u),lc),  is  the  associated  amplitude  of  the  wave.  As  t  -*  <»,  Eq.  (1)  may 
take  on  one  of  the  two  forms  Illustrated  in  Figure  1.  The  disturbance 
may  grow  in  time  at  every  point  in  space;  or,  if  the  speed  of  the  distur¬ 
bance  is  sufficiently  large,  it  may  propagate  away  while  growing  with  time 
and  leave  the  spatial  region  behind  unperturbed.  The  former  situation  is 
referred  to  as  an  absolute  instability  and  the  latter  as  a  convective 
instability. 


6 


The  Importance  of  making  the  above  distinction  should  no«r  be  clear. 

If  a  system  Is  oonveetlvel,  unstable,  a  steady  state  can  exist,  and  hence 

a  real  frequency  u,  and  the  corresponding  complex  «a,enumber-,ector  S(w) 

solutions  of  the  dispersion  equation.  The  original  oscillatory  «ave  Is 

being  amplified  as  It  propagates  along  In  the  system.  This  effect  «lll 

be  totally  washed  out  If  an  absolute  Instability  Is  present.  In  the 

situation  Of  absolute  Instability,  the  response  grows  monotonlcally  with 

time,  at  any  point  In  space,  until  It  Is  ultimately  llm'lted  by  nonlinear 
effects. 


PULSE 

AMPLITUDE 


PULSE 

amplitude 


FIGURE  1  SKETCHES  ILLUSTRATING  ABSOLUTE  AND  CONVECTIVE  INSTABILITIES 
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IV  INSTABILITY  ANALYSIS — CONFORMAL-MAPPING  METHOD 


One  basic  analysis  that  explains  the  Instability  theory  and  defines 
the  various  types  of  Instabilities  In  an  unstable  system  Is  accomplished 
by  using  the  dispersion  equation  together  with  the  appropriate  Initial  and 
boundary  conditions,  and  performing  a  conformal  mapping  from  the  complex 
0)  plane  onto  the  complex  Ic  plane,  or  vice-versa.  Our  task  In  this  section 
Is  to  present  a  detailed  study  of  the  conformal  mapping.  During  the  pro¬ 
cess  of  the  mapping,  the  growth  rate  for  the  wave  can  also  be  obtained 
when  the  system  Is  unstable. 


For  simplicity,  we  shall  consider  only  one-dlmenslonal  problems,  The 
one-dlmenslonal  analysis  actually  covers  an  Infinitely  uniform  system  or 
a  transverse  eigenmode  of  an  Infinitely  long  transversely  bounded  system. 

For  elementary  perturbations  propagating  as  exp  j  (cut  -  kz) ,  the 
response  of  a  linear  system  can  be  expressed  by  the  following  Laplace- 
Fourler  Inverslng  Integral: 

r(z,t)  = 

4it^ 


dk 


S(u),k)  J  (cut  -  kz) 
D(u),k)  ® 


(2) 


where 

8(0), k)  =  Laplace-Fourler  transform  of  the  source  function  S(t,z) 
D(a),k)  =  Dispersion  relation  of  the  system 

L  =  Laplace  contour  In  the  complex  u)  plane. 

The'  Integration  with  respect  to  k  In  Eq.  (2)  Is  carried  out  along 
the  real  k  line  as  shown  In  Figure  2.  At  this  stage  we  assume  that  S(u),k) 
■  Is  an  entire  function  and  also  that  D(u),k)  =  0  has  no  branch  points  In 


the  complex  k  plane.  Thus,  the  Integrand  of  Eq.  (2)  has  only  pole  singu¬ 
larities  in  the  complex  k  plane.  Using  Cauchy's  residue  theorem,  we  per¬ 
form  a  contour  integration  in  the  complex  k  plane  by  closing  the  k  contour 
in  the  lower  (or  upper)  half-plane  for  the  response  z  >  0  (or  z  <  0) .  The 
appropriate  k  contours  are  Illustrated  in  Figure  2  and  the  response  is 
given  by 


r^(z,a))  = 


z 

1=1,2, 


2) 


z  >  0 


(3a) 


n 


r^(z,co) 


J 


Z 

1=1  9. 


J^ojt  -  k  z^ 


m 


z  <  0 


(3b) 


where  k  and  k  represent  the  distinct  simple  zeros  of  D(u),k)  =  0  in  the 
n  m 

lower  half-plane  and  the  upper  half-plane,  respectively. 


The  above  choice  of  k  contours  for  z  >  0  and  z  <  0  is  subject  to  the 
radiation  conditions  at  |z|  -*  »  (boundary  conditions  at  infinity),  and  for 
t  >  0,  Eq.  (3)  represents  outgoing  waves  propagating  away  from  the  source. 


To  determine  the  Laplace  contour  for  integration  with  respect  to  o), 
we  Impose  the  initial  condition  that  response  is  zero  for  t  <  0.  Thus 
the  Laplace  contour  should  be  chosen  along  a  line  Im  cu  =  =  -a(a  >  0) 

below  the  lowest  pole  of  r(z,u)).  The  contour  should  then  be  closed  in 
the  upper  half -plane  for  t  >  0  and  in  the  lower  half -plane  for  t  <  0,  as 
Illustrated  in  Figure  3.  Since  r(z,u))  has  no  poles  in  the  region  below 
the  Laplace  contour,  there  is  no  response  for  t  <  0,  as  demanded  by 
causality.  For  t  >  0,  it  is  clear  from  Eq.  (2)  and  Figure  3  that  inte¬ 
gration  along  the  large  semicircle  yields  no  contribution.  Therefore,  the 
Integral  with  respect  to  u)  can  be  evaluated  by  summing  the  residues  at 
the  poles. 


9 


It  Is  now  clear  that  when  r(z,a))  has  poles  within  the  strip  0  2u>^2-a  , 
an  Instability  In  some  sense  may  exist  In  the  system.  To  Investigate  this 
Instability,  we  use  the  dispersion  equation  D(u>,k)  =  0  and  map  the  strip 
0  2u>^2-a  onto  the  complex  k  plane.  From  the  types  of  loci  for  the  corre¬ 
sponding  poles  In  the  complex  k  plane  we  can  determine  the  nature  of  the 
Instability. 

Generally,  one  does  not  wish  to  study  the  response  of  the  system  for 
all  time  (0:St<<»).  It  will  suffice  to  know  the  behavior  of  r(z,t)  as 
t  -*  00.  In  the  limit  of  t  -*  oa,  the  response  will  ultimately  be  dominated 
by  the  fastest-growing  component;  hence,  our  Interest  should  be  focused 
on  the  pole  closest  to  the  line  =  -a* 

A.  Conformal  Mapping  of  the  o)  Plane  onto  the  k  Plane 

Consider  the  situation  shown  In  Figure  4.  Let  u>su>  +Ju>  (u>  <0) 

o  or  ol  ol 

be  a  pole  below  the  real  u)  line.  The  corresponding  poles  In  the  k  plane 

are  determined  from  D(u>,k)  =  0,  as  shown  In  Figure  4(b).  We  now  let  u> 

ol 

be  less  negative — l.e.,  move  the  Laplace  contour  toward  the  real  u>  line, 
as  shown  In  Figure  4(a).  In  Figure  4(b),  we  Indicate  how  the  corresponding 
poles  In  the  k  plane  move  as  u>  ^  Increases  from  -a  toward  zero.  For  the 
case  of  no  double  k  root  of  D(u>,k)  =  0  considered  here,  two  types  of  loci 
are  possl.''.le:  The  loci  may  remain  In  their  half -planes,  or  they  may  cross 
the  real  k  line.  We  remark  that  as  a  pole  crosses  the  real  k  line  there 
will  be  a  discontinuous  Jump  In  the  values  of  r^(z,u>)  and  r^(z,cjo).  This 
Is  caused  by  the  transfer  of  a  residue  to  the  other  half -plane.  The  con¬ 
cept  of  analytic  continuation  should  be  used  to  derive  the  correct  response — 
l.e..  If  we  wish  to  determine  the  response  by  Integrating  In  a  region  above 
the  Laplace  contour  (o)^  =  -o  line),  (say,  along  the  o)^  =  0  line),  we  should 
also  deform  the  real  k  contour  In  such  a  way  that  It  only  Includes  residues 
at  the  poles  that  were  In  a  given  half-plane  when  values  of  o)  on  the  Laplace 
contour  were  used.  Thus,  the  analytic  continuations  of  r^(z,a))  and  r  (z,a)) 
are  consequently  obtained  by  following  the  deformed  contour  Indicated  In 
Figure  4(b). 
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It  is  now  clear  how  to  define  convective  instability  and  amplifying 

waves.  Several  cases  for  a  wave  competent  with  frequency  co  =  co  givine 

or 

rise  to  amplifying  and/or  evanescent  waves  are  Illustrated  in  Figure  4(b). 
To  determine  the  maximum  growth  rate,  the  dispersion  equation  should  be 
solved  for  all  real  cu  to  find  the  maximum  value  of  Im  k. 

In  the  situation  discussed  above,  the  dispersion  equation  has  no 

double  roots  of  k.  If  in  the  process  of  moving  the  Laplace  contour  toward 

the  real  o)  line,  D(a),k)  has  a  double  root  of  k  at  (o)  ,k  ).  In  the  vicinity 

s  s 

of  this  root  we  expand  D(a),k)  =  0: 


where  the  various  derivatives  are  to  be  evaluated  at  (co  ,k  ). 

s  s 

Since  D(a)  ,k  )  =  OD/3k)  =  0  (the  condition  for  a  double  root  or  a 

S  S  B’ 

saddle  point  of  k),  Eq.  (4)  then  yields  two  roots  at 


k  =  k 

±  « 


1/2 


(5) 


It  is  seen  from  Eq.  (5)  that  in  the  limit  of  o)  -  co  ,  the  two  poles 
k^  merge  with  each  other  at  k^.  This  situation  is  illustrated  in  Figure  5. 

Arbitrarily  associating  one  root  in  the  upper-half  k  plane  and  the 
other  In  the  lower-half  k  plane,  Eqs.  (5)  and  (2)  yield  the  identical 
response  for  z  >  0  and  z  <  0: 


- — « 


13 


FIGURE  5  PINCHING  OF  THE  CONTOURS  OF  r^(z.w)  BY  TWO  POLES 


It  1.  »«.  fro,  Eq.  M)  or  Eq.  (6)  th.t  <0  .  a,  1.  .  braroh  point  In 
the  0)  plane  and  that  the  contribution  at  o)  =  dominates  the  response. 
Note  that  as  o)  -•  cu^  ,  Eqs.  (5)  and  (6)  not  only  give  rise  to  the  coales¬ 
cence  of  two  poles,  but  also  result  in  the  pinching  of  the  contour  in 
the  k  complex,  as  shown  in  Figure  5. 


If  this  is  the  singularity  with  the  fastest  growth  rate,  we  sub¬ 
stitute  Eq.  (6)  into  Eq.  (2)  and  perform  the  Laplace  inversion.  The 
asymptotic  response  of  the  system  is  given  by 


r(z,t) 


anj  m 


|z|) 

e 


for  all  z  . 


(7) 


Thus,  the  effect  of  the  pinching  due  to  the  occurrence  of  a  saddle 
point  in  the  k  complex  plane  makes  the  response  infinite  at  cu  =  co  (note 
that  <  0)  for  all  z.  This  type  of  instability,  involving  merging 
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poles  in  the  k  complex,  is  termed  absolute  instability.  At  any  point  in 

space,  the  response  grows  monotonically  with  time,  until  it  is  ultimately 
limited  by  nonlinear  effects. 

It  may  be  remarked  that  coalescence  of  poles  that  de  not  pinch  the 
-entour  gives  a  finite  response  (no  absolute  instabilities). 

From  Eq.  (4),  »e  nete  that  OD/to)^  ^  0  (since  is  not  a  double 
roet  of  D(<„,k)].  The  condition  OD/dk)^  =  0  then  results  in  the  vanish¬ 
ing  of  greup  veleclty  v^  =  (&„/3k)^  at  point  (c^.k^).  Thus,  in  the  sit¬ 
uation  of  absolute  instability,  the  energy  of  the  respense  does  not  propa¬ 
gate  but  increases  monetonlcally  with  time  mlthln  the  entire  system, 

“•  Mapping  of  the  k  Plane  ente  the  m  Plane-Free - - 

An  alternative  procedure  to  determine  the  nature  of  instability  is 
suggested  by  the  fact  that  a  saddle  point  at  k^  in  the  k  plane,  with  cr  =ro 

•nd  <0^^  <  0,  naps  a  branch  pelnt  in  the  lower-half  or  plane.  Thus,  one 

can  now  reverse  the  previous  process  of  analytic  continuation  (in  the  k 
plane)  by  mapping  the  lines  parallel  te  the  Im  k  axle  from  the  k  plane 

onto  the  or  plane.  Eq.  (4)  or  (5).  for  the  occurrence  of  absolute  Insta- 
bility,  can  be  inverted  in  the  form; 


00  Ri  00 

s 


(8) 


Equation  (8)  shows  that  one  of  these  contours  in  the  complex  ao  plane 
necessarily  traces  a  cusp  and  thereby  locates  the  desired  frequency  ao  for 
the  absolute  instability.  ® 


However,  we  note  that,  from  the  previous  saddle-point  analysis,  a 
saddle  point  that  does  not  pinch  the  contours  from  the  upper-  and  lower- 
half  planes  results  in  no  absolute  stability.  In  order  to  differentiate 
between  the  types  of  instability  by  this  reverse  mapping  process,  it  is 

important  to  know  the  additional  properties  of  a  cusp  to  determine  an 
absolute  instability. 
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From  the  analysis  in  the  previous  section,  we  note  that  a  saddle  point 

k  ,  which  pinches  the  contours  and  results  In  an  absolute  Instability,  Is 
s 


produced  by  the  merging  of  the  two  poles  from  the  two  opposite  half¬ 
planes  and  the  pinching  of  the  contours.  In  fact,  the  two  half  k  planes 
can  be  considered  as  two  Rlemann  sheets  In  the  k  space  when  we  trace  the 
loci  of  the  roots  of  D(a),k)  =  0  In  the  k  space,  and  the  pinched  contour 
can  be  regarded  as  a  branch  cut.  The  corresponding  situations  In  the  o) 

plane  are  as  follows:  (1)  the  saddle  point  at  k  maps  Into  a  branch  point 

s 

o)  and  In  the  vicinity  of  cd  the  contour  forms  a  cusp,  and  (2)  the  cusp 
s  s 

starts  and  terminates  on  two  different  Rlemann  sheets  of  the  cd  plane 
(corresponding  to  the  pinching  of  the  contours  In  the  k  plane). 

We  now  summarize  the  procedure  as  follows: 

(1)  We  take  the  dispersion  equation  and  map  the  real  k  line  Into 

the  o)  plane.  In  general,  a  discrete  number  of  modes  u)  (k)  Is 

n 

thus  obtained. 


(2)  If  Im  u)  (k)  ^  0  for  n  and  k,  the  system  Is  stable  and  supports 

n 

only  evanescent  waves. 

(3)  If  Im  CD  (k)  <  0  for  some  n  and  range  of  k— say,  k  <  k  ^  k  — 

n  '  A  ^  B 

the  system  Is  unstable.  To  determine  the  Instability  type,  we 

map  lines,  k  =  constant,  with  abscissas  k.  ^  k  <  k  Into  the  u> 
r  A  B 

plane  and  find  a  cusp  that: 

(a)  Is  centered  at  a  branch  point  u)  with  Im  u>  <0 

s  s 

(b)  Starts  and  terminates  on  two  different  Rlemann  sheets 
of  the  CD  plane. 

If  the  above  conditions  3a  and  3b  are  satisfied,  the  Instability  Is 
absolute;  otherwise  It  Is  convective. 

A  qualitative  diagram  for  these  contours  In  the  u>  plane  Is  shown  In 


Figure  6. 


V  INSTABILITY  ANALYSIS  IN  WARM-PLASMA  SYSTEMS 


In  the  previous  analysis  we  assumed  that  D((i),k)  =  0  has  no  branch 
points  (and  thus  no  branch  cuts)  in  the  complex  k  plane.  It  follows 
that  there  is  only  one  form  of  the  dispersion  equation  for  the  entire 
k  complex.  When  kinetic  theory  is  used  to  Investigate  the  problems 
of  wave  Interaction  with  warm  plasma,  the  dispersion  equation  can  have 
branch  points  In  the  complex  k  plane.  The  wave^’varm-plasma  system  can 
then  have  more  than  one  distinctly  different  dispersion  equation, 
depending  on  the  regions  In  the  k  plane.  For  example,  the  dispersion 
equations  for  the  propagation  of  longitudinal  waves  In  a  collisionless, 
hot,  Isotropic  plasma  are  given  by 


D^(k,«i)) 


(9) 


where  denotes  unperturbed  velocity  distribution  in  one  direction, 

and  denotes  the  principal  part  of  the  Integral  with  the  Integration 
along  the  Im  v  axis  (note  that  Im  v  =  v^  =  (oo/k)^). 

The  branch  cut  in  the  k  plane  associated  with  Eq.  (9)  Is  shown  In 
Figure  7(a)  with  D^(k,(u)  valid  above  it  and  D  valid  below  It.  To 
calculate  the  response  r(z,(u),  one  is  faced  with  the  problem  of  deciding 
how  to  carry  out  the  k  Integration  through  the  branch  cut.  In  view  of 
the  branch  cut  not  being  extended  to  k  =  0  (this  is  due  to  the  fact  that 
no  particles  have  v  >  c,  the  speed  of  light  In  vacuum),  one  possible 
Integration  contour  is  Illustrated  in  Figure  7(b).  Now,  as  far  as 
stability  analysis  is  concerned,  it  Is  clear  that  the  branch  line  lies 
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below  the  k  real  axis  for  k  (=Re  k)  >  0.  Thus,  all  of  the  considerations 

r 

studied  in  Section  IV-A,  "Conformal  Mapping  u)  onto  the  k  Plane”  for 

determining  the  instabilities  are  preserved  if  one  uses  D  ((u,k)  fe” 

k  >0  and  D  (u),k)  for  k  <0. 
r  -  r 


(•I  BRANCH  CUT  FOR  W  ON  THE 
LAPLACE  CONTOUR 

FIGURE  7  BRANCH  CUT  AND  MODIFIED  PATH  IN  THE  COMPLEX  k  PLANE 

A  useful  criterion  for  determining  a  system  being  unstable  has  been 
developed  by  Penrose,  and  is  known  as  the  Penrose  Criterion.  To  describe 
the  Penrose  Criterion,  we  consider  the  dispersion  equation  for  one¬ 
dimensional  longitudinal  waves  in  a  uniform  non-Maxwell i an  plasma: 

09 
i 

2 


9f  (v)/3v 
o 


U)/k 


dv 


(10) 


The  Penrose  Criterion  states  that  the  system  can  be  unstable  only 
if  Z(u),k)  can  be  real  and  positive.  To  prove  the  above  statement,  we 
assume  that  perturbations  vary  like  exp  J  {.nt  -kz).  it  was  shown  in 
Section  IV  that  for  a  system  to  be  unstable,  the  dispersion  equation 
must  have  a  solution  with  Im  uj  <  0  for  real  k. 


Let  k  be  real  in  Eq .  (10).  It  follows  that  instability  exists  only 
if  Z(u)/k)  takes  a  real  positive  value  somewhere  in  the  lower  half  of  the 
U)/k  plane.  Note  that  if  we  assume  that  the  waves  vary  like  exp  -j  (cut  -  kz), 
then  Z(u)/k)  takes  a  real  positive  value  somewhere  in  the  upper  half-plane, 
permitting  the  existence  of  instability. 

Equation  (10)  defines  a  mapping  from  the  v  plane  to  the  Z  plane. 

If  instability  is  to  exist  in  the  plasma  for  a  given  distribution  f  (v) , 
solutions  of  Eq.  (10)  must  exist  for  which  Z  is  real  and  nop-negati^e 
when  Im  V  <  0  (k  is  taken  to  be  positive).  It  is  then  clear  that  the 
condition  for  instability  is  Just  that  the  region  of  the  mapping  of  the 
lower  half  of  the  v  plane  upon  the  Z  plane  by  means  of  Eq.  (10)  must 
contain  some  portion  of  the  positive  real  axis. 

The  above  criterion  is  particularly  useful  for  the  determination  of 
the  marginal  instability  (Im  v  -•  ,0)  for  an  unstable  system. 


VI  PROCEDURE  FOR  ANALYZING  INSTABILITIES 


It  is  useful  to  sununarize  a  few  steps  in  the  procedure  for  analyzing 
instabilities  [assuming  that  the  waves  vary  like  exp  J(a)t-kz)]: 

(1)  Given  a  dispersion  equation,  the  first  step  is  to  look  for 
(U-complex  roots  for  k  real.  If  there  are  no  roots  with 

<  0,  the  system  is  stable  and  supports  only  undamped  or 
evanescent  waves.  Under  these  conditions,  it  is  legitimate 
to  solve  for  k-complex  roots  with  u)  real. 

(2)  If  there  are  roots  with  u)^  <  0  in  (1),  the  system  is  unstable. 
We  can  perform  one  of  the  two  mappings  to  determine  the  type 
of  instability:  (1)  The  dispersion  equation  can  be  used  to 
map  u)  into  the  k  plane.  Plots  such  as  Figure  8  may  be 
obtained  showing  the  saddle  points  and  poles  for  discrimina¬ 
tion  between  absolute  and  convective  instabilities.  (2)  The 
dispersion  equation  can  be  used  to  map  the  real  k  axis  into 
the  (u  complex.  Plots  such  as  shown  in  Figure  9  may  be  ob¬ 
tained  to  determine  the  absolute  instability. 

(3)  Conformal  mapping  is  usually  a  very  tedious  process.  It  is 

thus  useful  to  apply  some  fairly  crude  checks  to  the  dispersion 
equation.  One  method  is  to  consider  the  solutions  of  the 
dispersion  equation  for  Im  (u  -  and  0.  If  all  the  k  roots 

of  D(k,aj^  -  j»)  =  0  lie  on  one  of  the  half  k  planes  and  the 
k  roots  of  D(k,u)^  -  JO)  =  0  lie  on  the  same  or  opposite  half¬ 
plane,  there  is  no  possibility  of  the  poles  pinching  the  con¬ 
tours.  Thus  the  system  is  either  stable  or  convectively 

unstable.  On  the  other  hand,  when  D(k,u)  -  Job)  =  o  has  k  roots 

r 
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2  3  4  5  '  B 

(e»  CONVECTIVE  INSTABILITY: 


FIGURE  8  QUALITATIVE  ILLUSTRATIONS  FOR  MAPPING  THE  LOWER-HALF 
«  PLANE  ONTO  THE  k  PLANE  TO  LOCATE  THE  INSTABILITIES 
AND  DETERMINE  THE  GROWTH  RATES 


on  different  i'lalf -planes,  it  is  possible  for  the  poles  to 

pinch  the  contours  as  Im  u)  increases  from  Consequently, 

we  check  for  double  roots  of  k  of  D(k,(i))  =0.  If  a  double 

root  k  exists  for  some  (u  with  Im  (u  =  ou  <  0,  then  there 

s  s  si 

is  an  absolute  instability  at  (u)  ,k  ). 

s  s 
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VII  EFFECTS  OF  DRIVING  SOURCE 


In  the  previous  sections  we  have  conducted  the  instability  analysis 

under  the  assumption  that  the  transformed  source  function  S(u),k)  is  an 

entire  function.  Thus  the  basic  characteristics  of  the  response  are 

determined  by  the  singularities  of  the  dispersion  equation  alone.  We 

now  consider  the  response  of  an  unstable  system  to  a  continuing  source 

at  spatial  origin  (z  =  0)  by  taking  the  source  function  S(t,z)  =  6(z) 

exp  When  is  real,  it  can  be  regarded  as  a  Fourier  component 

of  a  periodic  excitation.  On  the  other  hand,  when  u)  =  cu  +  jcu  is 

e  er  ei 

complex,  it  can  be  regarded  as  an  external  excitation  of  increasing 
((Bei  <  0)  or  decreasing  >  0)  amplitude  according  to  (exp  -  u)  t) 

exp  jfu  t. 
er 

The  analysis  described  in  the  previous  sections  for  evaluating 
the  asymptotic  response  proceeds  as  before.  However,  in  evaluating 
the  Laplace  integral  of  Eq.  (2),  we  must  take  account  of  the  singu¬ 
larity  of  S(u),k)  =  [j(u)  -(!))]  (i.e.,  the  pole  at  u)  =  U)  ) .  It  is 

e  e 

clear  that  if  o)  lies  above  the  lowest  relevant  branch  point  (u)  ,k  ) 

®  s  s 

of  D(u),k)  =  0  (i.e.,  if  >  min  then  the  asymptotic  response 

will  still  be  given  by  Eq .  (7).  On  the  other  hand,  if  u)  lies  below 

e 

the  lowest  relevant  branch  point  (i.e.,  if  u)  <  min  u)  ),  then  the 

ei  si 

asymptotic  response  will  be  governed  by  the  applied  source,  which  can 
be  evaluated  by  the  use  of  Eqs.  (2)  and  (3): 

Ju)  t 

r(z,t)  =  r(z,w  )  e  ® 
e 

t  -  ®  (11) 
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where  r(z,u)^)  is  given  either  by  (3a)  or  (3b)  with  exp  Ju)t  suppressed, 

and  the  sununation  in  Eq.  (3)  is  now  over  the  inodes  k  (u)  )  or  k  (u)  )  . 

n  e  me* 

Note  that  r^(z,u)^)  or  r^(z,u)^),  depending  on  whether  z  >  0  or  z  <  0, 

should  be  analytically  continued  to  determine  whether  a  particular 

root  k(u)^)  appears  in  the  response  and,  hence,  from  the  sign  of  k  (u)  ) , 

i  e 

whether  it  grows  or  decays  in  space  away  from  the  source. 

Hitherto,  in  this  type  of  analysis,  u)  has  been  taken  as  real 

e  ’ 

corresponding  to  steady  excitation,  and  the  analysis  is  applicable 
only  for  systems  that  are  not  absolutely  unstable.  In  this  case  the 
argument  given  above  leads  to  criteria  for  deciding  on  which  side  of 
the  driving  point  the  various  waves  k(a)  )  appear,  and  hence,  from  the 
sign  of  k^,  whether  they  are  amplifying  or  damped  waves.  The  extension 

to  complex  is  needed  for  bounded  systems  and  is  not  considered  in  this 
report . 

In  the  above  context,  it  is  physically  clear  that  even  if  the  system 
is  absolutely  unstable,  provided  the  drive  Increases  more  rapidly  than 
the  fastest  growing  absolute  instability,  the  asymptotic  response  will 
be  governed  by  the  source  and  not  by  the  natural  response  of  the  system. 
On  the  other  hand,  if  the  system  is  not  absolutely  unstable,  the  asymp¬ 
totic  response  will  again  be  governed  by  the  source,  even  though  it  is 
a  decreasing  one,  provided  it  decays  less  rapidly  than  the  slowest- 
decaying  natural  response  of  the  system.  In  either  case  the  criteria 
usually  employed  for  real  o)^,  for  determining  whether  the  roots  k((D  ) 
appear  in  the  response  for  z  <  0  or  z  >  0,  may  still  be  employed  for 
complex  CD  . 


26 


VIII  CONCLUDING  REMARKS 


In  the  present  report  we  have  presented  an  elementary  but  detailed 
analysis  on  the  theory  of  wave  Instability.  Our  analysis  Is  based 
on  a  conformal  mapping  through  the  use  of  the  dispersion  equation. 

This  method  Is  by  no  means  the  only  one  to  analyze  the  Instabilities. 
Other  approaches  such  as  those  based  on  an  energy  principle  or  the 
semi-quantum  mechanical  approach  can  also  be  useful  In  examining  the 
Instabilities  of  a  wave  system.  These  two  approaches  are  particularly 
useful  to  Investigate  the  Instabilities  for  a  nonlinear  system.  Never¬ 
theless,  the  conformal-mapping  method  presented  here  Is  commonly  known 
to  most  scientists  and  engineers  and  Is  most  easily  appreciated. 

The  distinction  between  absolu^^'*  and  convective  Instabilities 
has  been  made  clear  by  now.  We  should  emphasize  the  Importance  of 
applying  Instability  criteria,  rather  than  blindly  solving  dispersion 
equations  for  an  arbitrary  choice  of  complex  k,  real  U),  or  complex  u), 
real  k. 


In  addition  to  the  two  variables  o)  and  k,  a  dispersion  equation 
usually  contains  a  number  of  parameters.  These  parameters,  based  on 
the  models  assumed,  are  the  mathematical  representation  of  certain 
physical  quantities  of  the  system— e.g.,  finite  boundaries,  weak 
Inhomogenelty ,  and  collisions.  Under  the  assumption  of  a  wave-system 
model,  the  first  ste^  Is  to  derive  the  dispersion  equation  D(u),k,Qr^)  =  0, 
where  represent  the  various  parameters. 


To  Investigate  the  Instabilities  of  the  system,  we  can  perform 
the  conformal  mappings  through  the  use  of  D(u),k,cir^)  =  0.  Since  the 
solutions  of  D(U),k,Qr^)  =  0  and  the  various  loci  In  the  mapping  process 
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depend  on  the  parameters  0^,  It  may  be  possible  to  control  or  adjust 
certain  physical  quantities  of  the  system  In  order  to  obtain  the 
desired  Instability  for  the  wave  growth.  This  point  Is  particularly 
Important  In  the  design  of  plasma-lnJectlon  experiments  for  proving 
certain  feasible  applications. 
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